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Recent measurements of the reaction ^H(d,p)'^H in metallic environments at very low energies 
performed by different experimental groups point to an enhanced electron screening effect. How- 
ever, the resulting screening energies differ strongly for divers host metals and different experiments. 
Here, we present new experimental results and investigations of interfering processes in the irradi- 
ated targets. These measurements inside metals set special challenges and pitfalls which make them 
and the data analysis particularly error-prone. There are multi-parameter collateral effects which 
are crucial for the correct interpretation of the observed experimental yields. They mainly originate 
from target surface contaminations due to residual gases in the vacuum as well as from inhomo- 
geneities and instabilities in the deuteron density distribution in the targets. In order to address 
these problems an improved differential analysis method beyond the standard procedures has been 
implemented. Profound scrutiny of the other experiments demonstrates that the observed unusual 
changes in the reaction yields are mainly due to deuteron density dynamics simulating the alleged 
screening energy values. The experimental results are compared with different theoretical models 
of the electron screening in metals. The Debye-Hiickel model that has been previously proposed 
to explain the influence of the electron screening on both nuclear reactions and radioactive decays 
could be clearly excluded. 

PACS numbers: 25.45.-z, 25.60.Pj, 26.20. +f, 23.90.+W 



I. INTRODUCTION 



The cross section for nuclear reactions between charged 
particles at low energies is mainly determined by the pen- 
etration probability through the Coulomb barrier, which 
results in a steep exponential decrease towards lower en- 
ergies. At sufficiently low energies, however, this decrease 
is slowed down due to screening the Coulomb barrier by 
the inevitable presence of surrounding electrons. The 
electron screening was originally taken into account for 
nuclear reactions precedingin dense astrophysical plas- 
mas in the interior of stars [1| where the nuclear reaction 
rates can be increased even by many orders of magnitude. 
For laboratory investigations of nuclear reactions at very 
low energies, this eflFect was theoretically described [2| 
and experimentally observed in different fusion reactions 
on gas targets, e.g. The corresponding enhancement 
of the nuclear cross section could be explained by the 
gain of electron binding energies between the initial dis- 
tant atoms and the final fused atom. This was attributed 
to the raise of the kinetic energy of colliding nuclei and 
called electron screening energy. For the first time, the 
electron screening effect resulting from much more im- 
portant for astrophysical applications free electrons was 
investigated in the d+d fusion reactions taking place in 
metallic environments 0, H, 0] ■ The experimentally de- 
termined screening energies for some heavier metals were 
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one order of magnitude larger than the gas target value 
[3| and larger by a factor of about four than the theoret- 
ical predictions These experimental results were also 
confirmed by other authors 0, [13, [HI, [13, H [ll- 

Meanwhile, the electron screening effect in the d+d fu- 
sion reactions has been studied for over 50 different met- 
als and some insulators [l^.[T3.[l3| allowing, in principle, 
for a systematic study of the target material dependence 
of the electron screening energy. Unfortunately, there are 
some discrepancies between experimental values obtained 
by different groups They probably arise from some 
experimental systematic uncertainties connected with ox- 
idation of the target surface or with a high mobility of 
the implanted deuterons under beam irradiation, which 
can lead to unstable deuteron density profiles within the 
target. Both effects play a crucial role for the experimen- 
tal determination of the screening energies mUM- The 
basic quantity received from the experiment is the nu- 
clear reaction yield which is given for a thick target as an 
integral over the range of the projectiles Y = [na] dx 
with the target nuclei density distribution n and the cross 
section a. So deviations in the observed yield have the 
two principal causes: changes in the deuteron density 
profile and modification of the cross section, probably by 
the screening effect, which are merged in the integrant 
product [nu]. Some standard experimental difficulties 
have been already discussed in our preceding paper 
where an especially adapted data acquisition and analysis 
method, allowing us to discern between changes in n and 
cr, has been presented in a systematic manner. Based on 
this method, we report here some new experimental re- 
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suits and estimate experimental uncertainties of previous 
experiments. We additionally compare data obtained by 
different groups and discuss systematic errors of applied 
experimental and analytical methods. 

Prom the theoretical point of view, the large number 
of experimental data corrected for the discussed experi- 
mental uncertainties enables a comparison with theoret- 
ical predictions. The first ab-initio quantum mechanical 
calculation of the screening energy in a crystal environ- 
ment has been recently performed using realistic wave 
functions However, the results are still unsatisfying 
because of the very high demand for computational power 
limiting the model accuracy. Thus, the self-consistent di- 
electric function theory developed previously [l^ will be 
used here for the calculation of the screening energy con- 
tributions coming not only from free electrons but also 
from bound electrons of reacting nuclei and host met- 
als. Additionally, the interaction with the crystal lattice 
will be included. The theoretical results will be extended 
for comparison with the last experimental studies of the 
electron scree ning in nuclear reactions between heavier 
nuclei [H, [13, [Ull and in radioactive decays [H, [H, [Hi . 
On the other hand, it has recently been suggested that 
the enhanced electron screening can be explained within 
the classical Debye-Hiickel model [l^. The idea has been 
supported by an observation of the predicted inverse pro- 
portionality of the experimental screening energies to the 
square root of the absolute target temperature jH, HB]- 
As a consequence one could expect a magnification of the 
a and decay rates when radioactive sources would be 
put in metals at cryogenic temperatures. Even though 
the Debye screening cannot be applied to strongly cou- 
pled electron plasmas, as metals at moderate tempera- 
tures are, the suggestion has found much public interest 
[13, [13, HI, [2§| ■ Thus, both experimental and theoretical 
aspects of the temperature effect of the electron screen- 
ing will be subject of a critical discussion clearly showing 
the inapplicability of the Debye-Hiickel model for these 
issues. 



II. EXPERIMENTAL SET-UP, DATA 
ACQUISITION AND ANALYSIS 

The experiments have been carried out at an accelera- 
tor optimized for low energy beams. Fig. [1] illustrates the 
principal set-up and the data acquisition system. The 
accelerator consists of a radio frequency ion source, an 
acceleration line powered by a highly stabiHzed 60 kV 
supply and subsequent electric quadrupoles for focusing 
and a magnetic dipole for beam analysis. The beam im- 
pinges onto a Faraday cup just inside the target chamber 
where beam adjustment can be done without disturbing 
the deuteron density in the targets. A horizontal mag- 
netic steerer is used to defiect the beam onto the target, 
such removing neutral particles and contaminations car- 
ried along by the beam. A cylinder box set to a negative 
potential surrounds the target in order to suppress sec- 



ondary electrons. The isolated target holder is connected 
to a current integrator. The targets were disks made from 
different pure metals becoming self-implanted deuterium 
targets under the beam irradiation. Four Si-detectors at 
the laboratory angles of 90°, 110°, 130° and 150° were 
used for the detection of all charged particles, p, t, "^He, 
of the reactions ^H(d,p)t and ^H(d,n)'^He. The detectors 
needed to be shielded from the backscattered deuterons 
in order to prevent a congestion of them and the data 
acquisition system. Therefore grounded Al-foils of thick- 
nesses from 120 — 150 |J.g/cm^ were placed in front of the 
detectors. The thickness is sufficient to block deuterons 
up to 60keV while all other ejectiles could pass. The 
detector voltage pulses travel through pre-amplifiers and 
spectroscopic amplifiers. The signals are digitized by four 
ADCs in an embedded VME system connected to a com- 
puter which automatically integrates the proton lines of 
the spectra in fixed time intervals [tHI and records the 
four differential counting numbers iV [&) and the charge 
q of the integrated beam current at the target in a file 
which then can be further processed. An example spec- 
trum is shown in fig [2 all ejectile lines are clearly identi- 
fiable. Due to the anisotropic angular distribution of the 
ejectiles of the d+d fusion reactions even at the lowest 
energies, a total counting number is calculated pro- 
viding the tabulated function (cj) which is the basic 
quantity for the further data analysis. 

Correspondingly, the experimental reaction yield is 
given by 



e dq 



(1) 



where the number of impacting projectiles is already sub- 
stituted by their charge, e is the detector efficiency and 
z the charge state of the projectile. On the other hand 
the yield can be calculated for an infinitely thick target 
(regarding the projectile range R) by 
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with the number density of the target nuclei n and the 
cross section a. Unlike other chemical compounds the 
small hydrogen atoms are not trapped in firm chemical 
bonds with metals. The hydrogen density is not bound 
to a fixed stoichiometric ratio and can and indeed does 
change under ion irradiation. Changes in the yield may 
now originate from both the deuteron density and the 
cross section and need to be discriminated. The den- 
sity is here a function of the target depth, the projectile 
energy, the implanted charge, the beam fiux and other 
target material dependent and environmental conditions. 
The tabulated function N (q) provided by our data ac- 
quisition system makes it possible to retain the differen- 
tiation in ([T]) and thereby gain information on the charge 
development of a depth averaged density n (q) . So as- 
suming depth homogeneity of the deuteron density in ^ 
the depth x can be substituted by the projectile energy 
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Figure 1: Experimental set-up 



E with the stopping power differential equation [3C 
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where cm and cd are the stopping power coefficients in 
the metal and in hydrogen, no the appendant hydrogen 
density. Applying this substitution one arrives at a mo- 
tivation and an expression for the reduced yield 0, 0, 
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Since both the cross section in the metallic environment 
and the deuteron density are unknown the yield need to 
be set in relation to a known gas target cross section. We 
therefore chose the parameterization from [3l| because it 
has the highest precision. It forms together with the low 
energy function (-v/E) of the stopping power ^ the in- 
tegral in the denominator on the right hand side. The 
gray printed expression is per se a constant. So if the re- 
duced yield is not constant it is based on deviations of the 
prescribed progression in the cross section or in the func- 
tional dependence of the stopping powers or changes in 
the density. It is a sensitive measure for such deviations 
but the distinction of the possible reasons is a matter of 
reasonable interpretation. Fig. [2] shows plots of the re- 
duced yield at two different energies. One can see long 
term changes in the individual measurements indicated 
by the straight lines. These are attributed to changes in 
the deuteron density profiles scattered by the counting 
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Figure 2: Analysis procedure at the example of Zirconium at 
lOkeV 



statistics, of course. In contrast, the large discontinuities 
of the reduced yield at the switching of the beam energy 
result from a modification of the cross section. This is 
taken into account by the enhancement factor F {E) in 
(0]). Since the absolute quantity of the deuteron den- 
sity is unknown for the practical analysis a normaHzed 
enhancement factor is defined 



y{E) _ F{E) 



(5) 



viEo) FiEo) 
with the normalization energy Eq which is chosen to be 
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Figure 3: Exemplary results for the enhancement factor 
Kiorm. Screening enhancement for Zr theoretically described 
by the curve with the single parameter Ue. 



yield can be used to calculate a deuteron density esti- 
mate by solving ^ for n (q) and supposing F = 1 [H, 
Eq. (10)]. Only for this purpose the stopping power coef- 
ficients are explicitly required. A corresponding density 
plot for an initial implantation in Al is shown in Fig. [TJ 
The numbers above the gray boxes in Fig. [2] are density 
estimates for that areas. 

This is in brief the basic experimental procedure as of 
0, @, [lB| • For the study of the electron screening effect 
two experimental campaigns were executed. Since the 
special physico-chemical properties of the hydrogen com- 
pounds and the beam induced chemical reactions at the 
target heavily infiuence the obtained results 0, [lB|; the 
second more extensive campaign needed to investigate 
these interfering effects [iB, Sec. 4] which are sketched in 
a concise survey in the next section IIIII 



25keV for the monitor measurements. The gray rect- 
angles indicate the points from which the experimental 
error for Fnorm is inferred. Thus not only errors from 
the counting statistics but also from long term changes 
of the density are included. Results obtained for different 
projectile energies are displayed in Fig. [31 Assuming elec- 
tron screening as the reason for the increase of i^norm and 
adopting Ue as a kinetic energy shift parameter called the 
screening energy in the cross section [3] of the yield one 
receives H, @, O 
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for the screening enhancement factor of thick target 
yields [t^. The factor 2 arises from the CM-Lab- 
transformation. So F is an enhancement factor for thick 
targets in analogy to the enhancement factor for thin 
targets from Q 
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using the S-factor parametrization of the cross section 
with the Sommerfeld parameter rj in the second line 
and applying an approximation in the third line, which 
demonstrates its qualitative behaviour as a roughly ex- 
ponential increase for decreasing energies. The corre- 
sponding curve in Fig. [3] obtained for a fitted value of 
Ue supports the screening hypothesis. Our data analy- 
sis procedure is thus independent of the absolute value 
of the deuteron densities inside the targets and the stop- 
ping power coefficients which otherwise would introduce 
errors of 10 — 20%. The functional dependency of the 
stopping powers on the energy ^/E has been repeatedly 
confirmed, see [321 and references therein. The reduced 



III. EXPERIMENTAL SPECIALTIES AND 
PITFALLS 

The investigation of nuclear reaction cross sections on 
deuterium in metals should be performed at the lowest 
possible energies. This means that the composition of the 
topmost atomic layers of the metaUic target is of crucial 
importance because of the quickly decreasing range of the 
beam ions, considerably below 1 |a.m. This exactly is un- 
usual for experimental nuclear physics. The usual set-ups 
in experimental nuclear physics are constructed in high 
vacuum technology. But here the contained water vapour 
from the surfaces of all materials leads under ion impact 
to a progressing oxidation of the target metal because 
of the stronger electron negativity of oxygen in compari- 
son to hydrogen. Hence, hydrogen is contained in metal 
oxides only in segregation at low and unstable densities. 
Consequently, the oxidation diminishes and eventually 
destroys the screening effect with the growing thickness 
of the metal oxide layer. Carbon hydrides contained in 
HV systems pose another problem leading to carbon lay- 
ers on the target as will be discussed below. In such a way 
generated alterations in the depth profile of the deuteron 
density distribution in the target is the singular domi- 
nating error source for the observed enhancement and 
the inferred screening energies. Our vacuum system is 
made of aluminium with elastomer gaskets pumped by 
turbo molecular pumps with auxiliary oil lubricated two 
stag e rotary vane pumps and LN2 cooled cryogenic traps 
[la . Fig. 1]. A residual gas analyzer (RGA) was used in 
order to monitor the composition of the residual gas in 
the vacuum. This is here merely a concise presentation; 
for a more extensive description see (l6| . 

In accordance to the literature about HV systems the 
main constituent of the residual gas is water. Water 
vapour is due to its extraordinarily high dipole moment 
very adhesive to solids and is hence chemisorbed to sur- 
faces. Now under the ion irradiation several processes 
are enabled. Via heating and phonon excitation at the 
surface the beam provides the activation energy for dis- 
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Figure 4: Scanning electron microscopic pictures of target 
surfaces. Left: Symptoms of embrittlement for Al. Right: 
Beginning layer formation for Ta in island growth mode. 



sociative chemisorption of the water molecule, i.e. the 
protons are splitted off and the remaining oxygen radical 
forms a chemical bond to the metal atoms. Essentially, 
the same happens by direct impact excitation of the wa- 
ter molecule by the ions. The hydrogen implantation 
into the metal causes aside from the usual surface deteri- 
oration a in depth destruction of the crystal integrity of 
the material known as embrittlement which always oc- 
curs if the hydrogen loading rate is too high and not 
proceeding in thermal equilibrium [s^l- Thus, the sur- 
face is fractalized and the oxidation can progress into 
the bulk of the metal quickly creating a thick metal ox- 
ide layer. Fig. |4] contains as an example for it a picture 
of the surface of an Al target which turned into a sponge 
like structure. The rate of the oxidation process depends 
on the concrete form of the mutual interaction potential 
between the water molecule and the surface atoms, es- 
tablishing a material dependency. The energy supply of 
the beam enables these processes even for the noble met- 
als. Albeit generally spoken, more reactive metals apt 
more to oxidation and embrittlement while for the lat- 
ter the structural difference between the metal and the 
metal hydride is more important. Aside from the overall 
beam heating the energy of the projectiles is also impor- 
tant because lower energy projectiles are more effective 
at the surface [s^l- The partial pressure of water in HV 
is so high that there are ample suppHes for the surface 
chemical reactions. The hit rate of water molecules with 
a sticking coefficient of almost one is in comparable or- 
ders of magnitude as usual beam currents of 10—100 |J.A. 
This impHes a dependence on the ion flux, too. There 
are two counteracting processes: Sputtering and ther- 
mal or ion stimulated desorption. The sputtering yield 
of the Hghtweight deuterons is far too low in order to 
keep the surface clean with the resulting sputtering rate. 
One would expect that an increased temperature of the 
surface would increase the desorption rate of the water 
molecules. If the activation energy barrier for dissociative 
chemisorption of water is positive an increased tempera- 
ture yet proliferates the oxidation[73|- Similar is valid for 
ion stimulated desorption/chemisorption. Such again de- 
pends on the interaction potential but usually oxidation 
prevails. Unless UHV systems equipped for entire bak- 
ing are used the oxidation cannot be avoided. A deuteron 
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Figure 5: Effects of different surface compositions on the in- 
ferred screening energy for Ta. Ta-A has a small C-excess, 
Ta-E has slight C-traces, Ta-C a thick C-layer, Ta-D a thick 
MOi-layer. 



irradiation of only 1 C is enough to produce a consider- 
able metal oxide layer, see [la . Fig. 6]. There is, how- 
ever, a process that is nonetheless able to prevent oxida- 
tion: large carbon hydride molecules, e.g. backstreaming 
from the forepumps, can be physisorbed at the surface, 
cracked up and the carbon atoms can react with the oxy- 
gen radicals to carbon monoxide keeping in that way the 
surface clean. Differently from water, carbon hydrides 
are physisorbed to surfaces. The strength of this weaker 
bond increases with growing molecular mass. The ratio 
of absorption and desorption under the ion irradiation 
has similar dependencies. An evidence for this chem- 
ical surface reaction is the detection of a considerable 
CO fraction by the RGA which was below the detection 
threshold without beam irradiation [H, Fig. 7]. These 
processes were thoroughly explored by the regulated in- 
fusion of decane with monitoring feedback as the main 
part of the second experimental campaign. The surface 
can only be kept clean if the fraction of water and carbon 
hydrides in the residual gas are in an equilibrium which 
is of course also dependent on prementioned parameters. 
If the fraction of carbon hydrides is too low the surface 
will oxidize. If it is too high a carbon layer will build 
up. Both is essentially irreversible. Fig. [5] shows some of 
the results of these experiments for Ta demonstrating the 
high spread in the inferred screening energies depending 
on the surface composition which were verifled by elec- 
tron dispersive X-ray micro analysis (EDX). 

In order to limit the layer formation the totally im- 
planted charge was reduced [iB, Sec. 4.2]. For the anal- 
ysis a more sophisticated expression for the yield in ([H 
[5| was used based on a model of the target with three 
stacked layers [iB, Sec. 4.3]: The top layer consisting of 
either metal oxide or carbon, a deuterized zone of the 
metal and the bulk of the metal containing essentially 
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no hydrogen. Each can have different thicknesses and 
relative deuterium contents. The results for Ue in FigH] 
were obtained with only the additional parameter for 
the thickness of the deuterated zone in the metal in en- 
ergy equivalent units of the stopping [l^ . Sec. 5]. The 
differences for Ta-A and Ta-E are already considerable 
though the thicknesses of the surface layers were small 
and just started forming. Fig. |4] shows the beginning 
of the formation of a carbon layer starting from islands 
which will eventually cover the whole surface in concor- 
dance with experiences from thin film technology [33 |. 
Ta-C has already a relatively thick carbon layer which 
strongly reduced the screening energy. Just as the metal 
oxide layer does in Ta-D. Those layers were just thick 
enough in order to be included in the model and infer 
their thickness. The thickness of the metal oxide layer 
is 0.09 VkcV, which conforms to about 7nm. The cor- 
responding screening energy would be 433 eV [l^, Sec. 
5, Table 2]. 15 nm are enou gh t o let the screening en- 
hancement completely vanish [la . Sec. 4.3]. The deduced 
deuteron density is hardly affected and still in the vicin- 
ity of the stoichiometric ratio as the example in Fig.[6l[a) 
shows [l^ . Sec. 6, Fig. 13.(e,f)]. Much thinner surface 
layers already reduce the inferred screening energy con- 
siderably. So the real value for the screening energy of Ta 
is possibly around 400 eV [H, Sec. 5, Table 2]. Anyhow, 
the screening energy values ranging from 210 — 460 eV 
give an imagination of the systematic error originating 
from the surface layer formation. Carbon can achieve 
high deuteron densities but it does not show the electron 
screening effect as Fig.[7]proves. Thin deuterated carbon 
layers can, however, simulate a screening enhancement 
as inhomogeneous density profiles can do [H, Sec. 4.3]. 
Though the former could be excluded in our experiments 
but is a theoretical possibility when thin deuterated car- 
bon layers form on targets containing few deuterium in 
segregation as below. 

As already said, the metal oxide contains only few deu- 
terium in segregation. Those low densities are unstable 
and change under different conditions. At the example of 
a Na target with a very thick metal oxide layer the devel- 
opment of the calculated deuteron density is illustrated 
in Fig. [6l^b) (also [11, Fig. 13.(b-d)]). The density esti- 
mates are calculated from the reduced yield as previously 
described. Before the monitor measurement at 25 keV a 
measurement at a low energy had been taken. Then the 
density quickly decreased at 25 keV. Thereafter a mea- 
surement at 12 keV were started. Now, the density very 
quickly increased reaching a higher level than at 25 keV. 
But the discontinuity at the beginning was in the oppo- 
site direction. The density for the sequencing monitor 
measurement started once again at a high density which 
quickly decreased. The discontinuity at the beginning 
was once again in the wrong direction. So there is def- 
initely no screening in contrast to the positive case of 
Fig. [21 The quick shifts in the densities after the change 
of the implantation energy going to a 'saturation' level 
originate from a shift of the deuteron distribution depth 
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counter examples to the screening enhancement as in Fig. [2] 
(a) A medium thick layer obliterates the screening enhance- 
ment discontinuities at high densities, here at the example 
of an oxide layer on Ta. The full line shows the progres- 
sion from a refinement of the statistics by recalculation with 
an increased stepsize. Targets featuring low hydrogen binding 
ability hence allowing only for low and unstable densities with 
quick profile shifts: (b) A thick metal oxide layer overtopping 
the ion range on a metallic Na disk, (c) Heating vanquishs 
the hydrogen metal bond, here beam heating of a 7 vim Ta 
foil. 
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profile in tlie metal oxide linked to the different ranges 
of the ions [B, Sec. 6, Fig. 14]. With our method of 
recording a yield function Y{q) over the implanted charge 
we can recognize those shifts and reject them. If, how- 
ever, only the total yields of the long time measurements 
are regarded as in the usually appHed standard method 
(like in the other experiments discussed in Sec. lV Ap their 
comparison would erroneously lead to a screening inter- 
pretation. 

The same problem arises when working with low im- 
plantation densities below the stoichiometric ratio even 
when the metal oxide layer is negligible. Except for insuf- 
ficient implantation the density remains low if the ther- 
mal energy of the deuterons is higher than their chem- 
ical binding energy to the metal so that they can float. 
This applies mainly to transition metals with low ability 
to bind hydrogen (groups 6A-8A, IB) or if the metals 
are heated. An example for the consequences of heat- 
ing is shown in Fig. [6l[c) for a Ta-foil of 7 ii.m which was 
heated by the beam power. One observes the same be- 
haviour and no real screening enhancement. The density 
returns to an equal saturation level if the surrounding 
conditions are the same, i.e. same beam energy, current, 
target heat flow etc. Tests with a Au-foil showed an alike 
behaviour. The most effective heat transportation mech- 
anism in solids is the free electron gas. Cooling the target 
holder has little effect since the thermal resistance at the 
connection is very high. Besides from heating the density 
proflle of the deuterons in target materials with low bind- 
ing ability for deuterons (metal oxides, metals with low 
affinity to hydrogen, metals at high temperatures) is also 
changed by direct projectile hits and close phonon gen- 
eration at the target deuterons depending on the beam 
energy. Furthermore, the metal oxide as a thermal insu- 
lator will be considerably heated by the beam power. It 
is therefore preferable to use thick target disks at moder- 
ate temperatures with high densities. On the other side, 
cooling a target to very deep temperatures would trans- 
form it into a cryogenic trap accumulating water in thick 
layers on its surface prior to irradiation promoting the 
oxidation. The detailed investigation is covered in 

Summarizing, our data analysis method is independent 
of the absolute deuteron density and allows for the dis- 
crimination between changes in the reaction yield due to 
the density dynamics as in Fig. [6Ub,c) which are rejected 
and actual changes in the cross section which become 
manifest in the discontinuities at the edges of the mea- 
surements Hke in Fig. [21 That the discontinuities signify 
cross section modifications is further ensured by analyz- 
ing measurements which are taken in proximity of the 
stoichiometric ratio only, where changes of parameters 
like beam fiux and temperature have marginal infiuence 
on the overall deuteron density, at most. The error of 
-Fnorm is a couvolutiou of the error from the counting 
statistics and long term changes of the density. The use 
of high vacuum systems will inevitably cause the build-up 
of contamination layers. Thanks to our analysis method 
those layers can only diminish the inferred screening en- 



Table I: Screening energies 



Metal 




Ue in eV 


Ta 


0.9 


322 ± 15 


Zr 


2.1 


297 ± 8 


Al 


0.8 


190 ± 15 


Sr 


1.0 


350 - 800 


Li 


0.03 


^ 150 


Na 


0.03 


c 


Pd 


0.3*' 


313 ± 2 


C 


d 






"Approximate average deuterium contents in relation to the num- 
ber density of the metal 
'The initial implantation was deliberately prematurely aborted. 
"^An oxidation layer impeded the determination of Ue 
''Carbon density unknown. See text. 



ergy since feigned enhancements due to density dynam- 
ics get rejected [iB, Sec. 6]. The utiHzation of carbon 
hydrides embanks the layer formation enabling the re- 
sults in HV at all [l^ . Sec. 4.2]. Indeed this is a difficult 
and labile equifibrium in the residual gas. So layers are 
present, which were examined by EDX allowing for a rel- 
ative measurement of element abundances [H]. But the 
thickness can hardly be quantified because of the fractal 
structure of the target surfaces (e.g. Fig.[4l). Though the 
model suggests that 15 nm are sufficient to completely 
dispose of the screening enhancement. All in all, the 
obtained screening energies represent lower limits to the 
real value. The magnitude of the dominating systematic 
error from the unknown layer thickness can be assessed 
by the measurements in Fig. [H 



IV. RESULTS 

A. Experiment 

The experimentally determined results for the screen- 
ing energies are summarized in Table [B The values from 
the campaign I are in the upper part of the table. In the 
lower part of the table are accessory results from the cam- 
paign II. In the second column of the table are the ratios 
of the deuterium number density to that of the hostmet- 
als. Since the deuteron density can and does vary during 
a measurement these values are estimated averages. The 
values for Strontium and Lithium are heavily impaired 
by the layer formation due to the high reactivity of both 
metals, for Li more than for Sr. Such expressed itself 
as strong variations in the deuteron densites and accord- 
ingly in the reduced yields during the course of the mea- 
surement, leading to ambiguous values for the disconti- 
nuities of the reduced yields. So these screening energies 
should be regarded as estimations, at best. The results 
were obtained utilizing the equilibrium in the residual 
gas in order to keep the target surface clean which was 
subsequently verified by EDX analysis (Sec. [Till and more 
detailed [11, Sec. 4.2]). 
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U = (313+2) eV 

^ °= (1.461+0.014) ^IkeV 
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with about 5-10^^cni~^. Howbeit, these are only esti- 
mates since the carbon densities are not known and as 
a result the correct stopping power coefficient which is 
required ^ (igI . Eq. (10)] neither. Anyway, the resulting 
enhancement factors show no significant disagreement. 
Thus, carbon films present no signs of electron screen- 
ing. These results are Hsted in the lower part of Table [B 
The highly reactive metal natrium corroded so easily 
that only low deuteron densities could be achieved and 
no screening was visible. Two tests with Y and Er led to 
thick metal oxide layers, too. Different to the other met- 
als the concomitant analysis of the '^He spectral peak re- 
vealed in both experimental campaigns for Li, Na and Sr 
a significant suppression of the neutron reaction channel 
and a simultaneous alteration of the angular anisotropy 



1.3- 
1.2 - 
1.1 - 
1.0- 
0.9- 
0.8- 



10 15 20 25 30 35 40 45 
Er(keV) 



Figure 7: Measured values of Fno 
different compositions 



for Pd and Carbon in two 



The first plot in Fig. [7| is a measurement on palladium 
with roughly equal residual gas conditions as for the Ta- 
measurements |l6l . Sec. 5]. The totally implanted charge 
was limited for the same reason, i.e. layer formation, as in 
the Ta-measurements of Fig. [5l The beam spot contains 
traces of carbon specifically some dark stains [iB, Fig. 
8]. The other two plots are the experimental prove that 
carbon has no screening enhancement. The drop to the 
lower energies originated from a lower deuterium content 
in the upper layers of the targets. This drop can also 
be caused by the voltage drop in the plasma inside the 
RF-ion source [IHl [iB, Sec. 2] which has a higher impact 
for lower energies relative to the monitor measurement 
at 25 keV. The two carbon targets were preparated with 
different methods. The first one was made by deposition 
of soot from a ethine fiame on a backing plate. A fiame 
of ethine (C2H2) burning with insufficient oxygene sup- 
ply produces very pure carbon. However, the material is 
amorphous and rather fiuffy. Accordingly, the deuteron 
density reaches only values of about 1.5 • 10^^ cm~^. The 
second target is a carbon film produced by the irradiation 
of aluminium with high decane pressure. That in such 
way deposited carbon was compactified by the impacting 
beam ions while forcing it to adopt the lattice structure of 
the substrate to a certain extend [s^. Hence the density 
of the carbon atoms is higher, so is the deuteron density 



B. Theory 

From the theoretical point of view the deuterized met- 
als can be treated as a strongly coupled plasma [§]. 
Since the velocity of reacting nuclei is smaller than the 
Fermi velocity, the electron screening effect corresponds 
to a static polarization of surrounding conduction and 
bound electrons. Consequently, the electrostatic poten- 
tial energy between reacting nuclei of charges Zi and Z2 
shielded in a metaUic medium can be described within 
the self-consistent dielectric function theory [l^ : 



Vir) 



Z^Z^e^ f 4Tre(pi (q) eip2 (q) 



{2TTf J e^{q)ec{q)q^ 

r-,0 Z\Z2e^ 

> 'J-pol 



exp (iqr) dPq 
(8) 



The wave-number dependent dielectric functions ei, and 
Ec describe polarization of valence and core electrons of 
host atoms induced by a charged impurity taking into 
account the short range electron correlation and the 
exchange interaction between electrons (for details see 
<I> (r) and tpi (q) functions are the screening func- 
tion and electronic charge-formfactors of reacting nuclei, 
respectively. At small distances (applicable for nuclear 
reactions and decays) the potential energy can be ap- 
proximated using the energy independent polarization 
screening energy Up which scales with the product of the 
charges of the involved nuclei. For the d+d reactions we 
used the self-consistent charge form-factor ip (q) within 
the Thomas- Fermi approximation [H, [Hi : 



ip{q) ^1- z- 



zq- 



Here, the Thomas-Fermi wave number k'^p 



(9) 



Qire'^n/Ep 



has been applied; n and Ep are the electron number 
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Figure 8: (Color online) Comparison between experimental 
and theoretical screening energies 



density and the Fermi energy, respectively. The num- 
ber z corresponds to the fraction of electrons bound to 
deuterons and is for metals close to unity. Since we are in- 
terested in the evaluation of the strongest possible screen- 
ing effect, we uniformly set z = 1 for all target materials. 
In the absence of screening = Ec = 1 and z = 0, V (r) 
reduces to the bare Coulomb potential (<& (r) = 1). 

In the metaUic lattice, besides electrons also positive 
ions can contribute to the screening of the Coulomb bar- 
rier between reacting nuclei. This effect, called cohesion 
screening can be calculated as a gain of the poten- 
tial energy of two deuterons in the lattice field of the host 
metal compared to that of the helium atom produced in 
the fusion reaction. To calculate the potential energies 
we used the universal ion-ion interaction given by Ziegler, 
Biersack and Littmark [3§|. For a rough estimation of the 
cohesion screening energy Ucoh, we calculated the poten- 
tial energy gain resulting from the surrounding 12 host 
atoms assuming the same fee crystal structure for all tar- 
get materials investigated. The cohesion screening is a 
slowly increasing function of the atomic number. The 
total screening energy is the sum of both contributions 

Ue = Upol + Ucoh- 

The results of the theoretical calculations obtained for 
the d+d reactions taking place in different metallic tar- 
gets are presented in Fig. [8] together with our experi- 
mental values. The electron screening energies moder- 
ately increase with the atomic number of host atoms 
reaching for heavier nuclei the value of about 300 eV. 
The experimental target material dependence agrees with 
the theoretical expectations. However, the experimental 
screening energies are larger by a factor of about 2 com- 
pared to the theoretical values. Since the experimen- 
tal screening energies obtained for insulating materials 
are much smaller (< 50 eV) and taking into account 
that the screening contributions arising from polariza- 
tion of bound host electrons and cohesion should be sim- 



Figure 9: (Color online) Experimental and theoretical polar- 
ization screening energies versus the electron gas parameter 
rs- For comparison the Thomas- Fermi screening of the elec- 
tron gas is presented. 



ilar for both metaUic and insulating targets, we can con- 
clude that the enhanced screening effect results from con- 
ducting electrons. Thus, for a comparison between dif- 
ferent target materials the electron-gas parameter rs = 
[3/ (47m)]'' /ao where n and ao are the valence electron 
density and the Bohr radius, respectively, is much more 
suitable. Using this parameter, the experimental polar- 
ization screening energies obtained by subtraction of the 
theoretical cohesion contribution are displayed in Fig. [9] 
together with the theoretical polarization screening en- 
ergies. Now, the quality of the theoretical description is 
much better visible. In contrast to the simple Thomas- 
Fermi model [3^ providing for free electrons a smooth 
dependence of the screening energy given by Utf = 
ZiZae^ [4/ (^o)]^ {Sn^ny = 2ZiZ2e^ [9/ (4^^)] i 
the dielectric function theory describes fluctuations of 
the experimental polarization screening energy very well. 
The fluctuations result from the polarization of bound 
(core) electrons which contribution to the total screen- 
ing energy depends very strongly on their binding energy 
If the bound electron contribution is eliminated 
from the experimental polarization screening energies we 
get experimental values for the free electron polarization 
which can be parameterized by a smooth dependence on 

rs Upoifj = ZiZ2(250 ± 20) eV/r|. This result can be 
used for an estimation of the free electron contribution 
in metallic environment to the screening energy in reac- 
tions between nuclei with charges Zi and Z2. Different 
to the d+d reactions, the contribution coming from elec- 
trons bound by heavier reacting nuclei is much larger 
and should be included separately. This can be calcu- 
lated as the gain in electron binding energies between 
distant atoms and the final united atom. Similar results 
are to obtain using the Thomas-Fermi model leading to 
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Ue^b (TF) = l.lSZiZse^ (^Zf + zj) ' /oq In the 

case of heavier nuclei the cohesion screening energy can 
be neglected, since the strength of the interaction with 
the lattice atoms increases much weaker than the product 
Z1Z2. Thus, the total screening energy is only the sum of 
the free electron and bound electron contributions. The 
same estimation can also be applied for radioactive a and 
P decays 

The dielectric function theory does not predict any 
temperature dependence of the polarization screening en- 
ergy unless the electron density of the target material re- 
mains constant and the projectile velocity is smaller than 
the Fermi velocity. That is typical for a strongly coupled 
plasma. For velocities higher than the Fermi velocity the 
electrons are not able to follow the ions and the electron 
screening gets weaker. In this limit of a weakly cou- 
pled plasma (Debye-Hiickel limit) the screening length 
becomes larger than the mean atomic distance and clas- 
sic description of the electron screening is applicable. The 
screening energy is inversely proportional to the square 
root of the kinetic energy or equivalently of the plasma 
temperature {Ue ^ ^ '^)' analytical formula 
connecting both limits has been derived by Lifschitz and 
Arista [4l| for the stopping power of moving ions in the 
electron gas and can be applied for the electron screening 
in nuclear reactions [l3|. Thus, the velocity dependence 
of the screening energy can be given as follows: 
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Ta 



-strong screening 
■weak (Debye) screening 
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Figure 10: (Color online) Screening energy dependence on the 
projectile energy. The Debye screening is applicable only for 
deuteron energies larger than the Fermi energy (56 keV for 
Ta) or equivalently for plasma temperature larger than the 
Fermi temperature (1.8 ■ 10^ K for Ta). 



A. d-|-d Experiments 



u, 



ad 



■In 



V + Vp 



V — Vp 



(10) 



where Udyn and Uad denote dynamic and adiabatic 
screening energies, respectively. The Fermi velocity vp 
depends on the electron density and therefore is charac- 
teristic for the target material. The above relation calcu- 
lated for the d+d reactions in the Ta environment is pre- 
sented in Fig.[TOl Additionally, the energy dependence of 
the Debye-Hiickel screening is shown. It is visible that the 
electron screening can be described by the Debye-Hiickel 
theory only for projectile energies higher than the Fermi 
energy (the Fermi energy of deuterons in Ta amounts 
to about 56 keV) or equivalently for temperatures higher 
than the Fermi temperature (for Ta ^ 1.8 • 10^ K). Thus, 
in the cases discussed here, the Debye-Hiickel screening is 
not applicable for both nuclear reactions and radioactive 
decays. 



COMPARISON WITH OTHER 
EXPERIMENTS 



In view of the augmented information provided by our 
differential analysis method and experimental procedure 
the results of other groups will be discussed. 



In Fig.[TT]an overview of screening energy results and 
appendant deuteron densities from other experiments is 
plotted. All were carried out in high vacuum systems 
hence suffering from the same progressive oxidation pro- 
cess under ion irradiation with the inherent problems as 
of Sec, mil and [l^ . A quick glance already shows that 
the screening energy results are pretty much scattered 
not revealing a pattern. But in conjunction with the 
deuterium metal ratio aka the deuteron density peculiar- 
ities become evident. Our high screening energy results 
(Table [H were achieved at high absolute densities in the 
proximity of the chemical stoichiometric ratio where the 
ion beam flux has no influence on the target deuteron dis- 
tribution. Whereas the high screening results of the other 
groups were exclusively attained at low deuteron densi- 
ties 10~^ -~^ below the metal number density. Comple- 
mentary high densities did not yield enhanced screening 
in those experiments. Both classes of screening results 
are associated to groups in the periodic table exposing 
the chemical relationship with respect to the surface re- 
actions and hydrogen binding ability of the targets as de- 
scribed in Sec. mil Particularly manifest is that in Fig.fTTl 
for the group 3A metals including the lanthanoides which 
have low screening values at high densities and conform 
to the counter example case of Fig. [6l[a). But also the 
transition metals show three clusters of high screening 
results at low densities in Fig. [iT] corresponding to the 
case of Fig. ^h) and Fig. [6l[c) respectively. Such will be 
substantiated in the following. 
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Bochum 2002 
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Figure 11: Overview of screening experiment results. Legend: Garching 1990 Tohoku 1998 2002 [T3|; Bochum 2002 
[ia |. 2003 2004 [l3|- Bottom: Screening energies Ue- Top: Deuterium to metal ratio x. The values for x of were 
estimated from Fig. 2 therein. The values of |l4l| are the data base; data points from [13, are included if they differ, only. 



The Garching Experiment 



results. 



The first accelerator experiment with the aim to search 
for modifications of the cross section in the d+d fusion 
reactions caused by the metallic environment was done 
on Ti [i^l- No enhancement could be observed. The 
measurements were performed on a 3 \mi Ti foil fixed 
in a copper target holder frame with flow channels for 
LN2 cooHng and a thermocouple for temperature deter- 
mination. No effort was made to specify the deuteron 
density in the target. Instead a fixed value from mate- 
rial research was adopted which is inadequate. Due to 
the deep cooHng of the target water is accumulated on 
its surface, which produces under ion irradiation a con- 
siderable oxide layer In addition, the beam current 
up to 0.1 mA leads to a distinct temperature increase in- 
side the de-acceleration volume of the ions in the thin 
foil which will also alter the density profile away from 
the supposed unit value. All further measurements on Ti 
in Fig. [11] resulted in very low screening values with den- 
sities in proximity of the chemical stoichiometric ratio. 
The higher the deuteron density, the lower the screening 
value. Ti is chemically very similar to Zr, both belong 
to the group 4A. From our experience Zr oxidizes read- 
ily. So a relativelythick metal oxide layer complying to 
the third case in |l6l . Fig. 13.(e,f)], Fig.[6lja) explains the 



2. The Tohoku Experiments 

The results of [1, are based on the analysis of the 
total yield of the proton measurements Eq. (1)] with 
one detector at the lab-angle 90° and at the projectile 
energy Ed S [2.5, 10] keV 



Yt (Ed) = sN^ 



Ed 1 



dE 



(11) 



after depth energy substitution (Sec. [ill with the proton 
detection efficiency e a cross section parameterization a 
of ji^l and finally the target deuteron number density 
A^D which is presumed to be constant for all energies 
and ranges. With the stopping power relation the ad- 
ditional error of the stopping power coefficients is intro- 
duced. In order to determine and observe the density 
value repeated monitor measurements were performed at 
lOkeV. The density was then calculated from the yield 
Yt (lOkeV) by solving ifTTj) for A^d with the supposition 
that the screening enhancement is there negligible. Ac- 
cording to the not unambigious text ([l3| and suitable 
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back references) for the quantification of the enhance- 
ment and extraction of the screening energy the yields 
are normalized to the experimental one at 10 keV thus 
becoming independent of the actual value of Nu : 



(E) 



fjE) -YtjE) 
Yt (lOkeV) 



Y{E) 

iaare {E) 



cxp m-j {Ecu) 



_Ue 

E, 



CM 



(12) 



(13) 



where the theoretical expression for the thick target en- 
hancement factor / is simply adopted from the approxi- 
mated term in ([7]) of [2| for the enhancement of the cross 
sections in thin targets. This expression is, however, de- 
rived for cross sections based on an increase of the ef- 
fective projectile energy. It might be introduced in the 
integrant in (fTTj) at the most and must not be pulled 
out of the integral in that manner. A thick target en- 
hancement factor should retain the energy integration 
like in The approximation as in ^ is only valid 

for -EcM ^ Ue which is not longer fulfilled by the given 
experimental circumstances with beam energies of some 
keVand screening energies of several hundred eV. This 
also means that the supposition of a neghgible screening 
enhancement at lOkeV is not valid either. Moreover, the 
term lfT3|) diverges for energies approaching zero. While 
the invalid approximation leads to an underestimation of 
the derived screening energy the neglect of the enhance- 
ment at 10 keV effectuates a gross exaggeration because 
the curvature of the enhancement curve must be greater 
in order to describe the steeper slope of the data (analog 
to the difference in the curves #1 and #2 in [JJ, Fig. 
12]). The inclusion of measurements taken at higher en- 
ergies would have revealed such. The deuteron density 
value obtained at lOkeV ought be heavily altered, too. 
The target holder was cooled with LN2. The constancy 
(and thence concluded the homogeneity) of the density in 
the target was investigated by measurements at lOkeV 
with target heating by different beam flux or a mounted 
heater in the case of Pd in the interval [170, 2301 K deter- 
mined with a thermocouple. The results in |lOl . Fig. 1] 
show a strong dependency of the 'saturation' density on 
the target temperature and material with a considerable 
general decrease with raising temperatures. The density 
descends from Ti over Au, Fe, Pd to PdO. Conspicuous 
are the differences in the deuteron densities between [§| 
and [l3| for Ti, Au and Pd which are almost one or- 
der of magnitude while the corresponding screening val- 
ues accord within their errors (Fig. [TT|) though the latter 
Ue are generally higher. This discrepancy remains unex- 
plained in [l3|- While Au and Fe do not build up firm 
bonds to hydrogen the achieved densities are proportion- 
ally higher at these deep temperatures. A deliberately 
produced 30 nm thick PdO layer on a Pd target in [10| 
yielded an especially high screening energy with an espe- 
cially low density. Such a thick PdO layer would show 
quick shifts in the deuteron density profiles with higher 



averaged densities at lower projectile energies like in the 
second case in [H, Fig. 13.(b-d)], Fig.[6l[b) when chang- 
ing the projectile energy and using the differential anal- 
ysis method. So this large screening is simulated by the 
density alteration during the total yield measurement. 
Whereas the density in the Pd target of [l3| is notice- 
ably low which points on the formation of an oxide layer 
distinctly thicker than on the Ti target approaching the 
PdO target. This means that the screening value is also 
rather generated by density dynamics and the agreement 
with our value is by chance. Albeit our Pd value was ob- 
tained from measurements with a limited total ion dose 
and still growing densities prior to saturation in order 
to minimize surface contaminations the highest density 
of all experiments was achieved. While the targets as 
described in [l3| are thick enough ~ 1 mm to guarantee 
an effective heat transport in the bulk of the material by 
the electron gas, the heterogeneous target Au/Pd/PdO 
with a total thickness of only 60 iJ.m (thereof 0.1 (a.m Au) 
is too thin therefore leading to a considerable tem- 
perature increase in the beam stopping volume which is 
to this extent not detectable by any outside mounted 
thermocouple. So, the observed high screening energy of 
(602 ± 23) eV can be explained by the shifts in the den- 
sity profile due to elevated temperatures like in Fig. [6^0), 
[ig . Fig. 13. a] and the heterogeneity of the target and 
accordingly the density. 

In order to explain the ascertained relation between the 
screening energy and the density depicted in [l3. Fig. 4], 
i.e. high screening comes along with low deuteron densi- 
ties, the concept of a deuteron 'fiuidity' was introduced in 
[lo| where fiuid deuterons and conduction electrons are 
to behave Hke a hot plasma. But in palladium oxide there 
are no conduction electrons. In view of the stated density 
dynamics this explanation is decrepit. The explanation 
by density dynamics is also sustained by the significantly 
larger standard deviations of the repeated density mea- 
surements at 10 keV for targets with low densities in 
Fig. 2]. Indeed, the saturation density in our experiments 
returns to the same level for the same conditions but with 
higher deviations. In order to avoid the observed temper- 
ature changes of the deuteron densities in the targets the 
beam current was adapted in such a way that the power 
input into the target was kept constant. Although in [l3| 
is admitted that this procedure does not keep constant 
the power density due to the stopping power relation the 
mobility of the deuterons is not only infiuenced by the 
indirect ambient temperature but also by direct ion in- 
teraction and changes in the distribution of the stopped 
projectile deuterons. The authors conceded that they 
were not able do detect possible short time changes in 
the proton counting rate. With our differential method 
we did not observe any discontinuity belonging to screen- 
ing on oxidized targets with low absolute densities being 
independent of the actual beam current and power in- 
put. The trustiest screening value seems to be the result 
for Au in [§| obtained at a high density which conforms 
to our test with a Au foil at a very low density yielding 
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no enhanced screening. But under ion irradiation even 
noble metals can oxidize. 

Postscript: Very recently the experiments have been 
continued using the same accelerator set-up, procedure, 
and analyzing method as above [i^l; consequently the 
same considerations apply. As before the target holder 
was cooled with LN2 further approaching the boiling tem- 
perature of nitrogen. For the Sm target a screening en- 
ergy of (520 ± 56) eV was deduced [79||. But not even an 
estimate for the deuteron density was given unHke in the 
previous puplications. The high statistical errors in [3, 
Fig. 5] and the vague statements regarding the deuteron 
'fluidity' obtrude the inference of a low density with the 
same consequences, too. 

3. The Bochum Experiments 

The largest data set of screening energies is provided 
by [H, flaTfl^l . The applied experimental procedure and 
data analysis method explained in [ill, is exactly the 
standard strategy in nuclear astrophysics as described 
in [4^ including measuring relative excitation functions 
with normalization to known cross sections. As such it 
is a step back behind [§| where already concessions to 
the special situation of hydrogen in metals were made. 
Again just the total yield for the thick target of the mea- 
surement of the protons with 4 detectors at the polar 
lab-angle of 6* = 130° at the energy G [5, 30] keV 
was determined. The total yield Y {Ed,0) was repeat- 
edly taken at flxed energies with equal stepsizes A of 
0.5 keV and 1.0 keV for > lOkeV. Thereof an energy 
differentiated yield Y'{Ed,0) is calculated lfT4|) in order 
to extract the cross section JTS]) [ll|, Eq. (5), (7), (8)][83|: 

Y'{Ed,9) - iY{Ed,9)^Y{Ed-A,d))/A (14) 
= a{Eee)eesiEdr^ x (15) 
nKa{Ed,0)W{Ed,e) 

a^const. 

The energy integration prior to ifTS]) vanished by means 
of the mean value theorem of calculus leaving behind 
the integrand to be evaluated at the effective energy 
EcS G {Ed ~ A, Ed) where one half of the yield is at- 
tained. Except for a for all other factors in ifTS]) is as- 
sumed that their change within the energy interval can 
be neglected. Moreover, the angular terms are collected 
in the factor a which is supposed to be constant for the 
whole energy range of the measurement, fl is the solid 
angle of the detectors, Kq its transformation to the CM- 
system and W the angular distribution of the reaction 
yield, e is there the stopping cross section, i.e. the en- 
ergy loss per particle areal density eV/ (atoms/cm^) and 
not the Hnear stopping power ^ kcV/|J.m. The effec- 
tive stopping cross section is assembled from the one for 
deuterium and the host metal [ill . Eq. (9)]: 

e^s{Ed) = eu{Ed) + xeuiEd) (16) 



with the metal atom fraction M^jD. Thus, the depen- 
dence of the composition of the target is completely shuf- 
fled from the stopping factor into x. Consequently, the 
deuteron density described by x is forced to be flxed for 
all projectile energies, the range in the target and the 
whole measurement series on the target. For the deter- 
mination of the absolute value of the cross section in ifTS]) 
X was scaled to a known cross section for gaseous deu- 
terium [3| at = 30 kcV. This means that the deduced 
uniform deuteron density for the whole measurement se- 
ries is only dependent on the one value at 30kcV and 
only there validated, at most. So this method is even 
less sensitive to changes of the density during the course 
of the experiment than [§|. Then the S-factor is calcu- 
lated. The screening energy is obtained from another 
flt to the S-factor data with three parameters of the ex- 
pression in the second line of ^ together with a linear 
S-factor function. Furthermore, additional error sources 
were introduced without need by sticking to the standard 
procedure: repeated yield differences at flxed energies, 
introduction of the effective energy, the stopping power 
coefficients, S-factor computation. The errors of the com- 
puted S-factors are said to be dominated by the spread in 
the yields Y' from various runs ( [ill . Sec. 4, p. 380], [l^ . 
Sec. 2, p. 195]), i.e. the yields were repeatedly measured 
with stepwise increasing and decreasing beam energies. 
(This implies that the errors from different F-values are 
distinctly higher than the corresponding statistical er- 
rors, which can be seen in [H, [ij. Fig. 1].) As can be 
seen from the position and errorbars of the datapoints 
m [13 ■ in . Fig. 1] the differences of the 1"- values must be 
significant. It complies with our experience that the den- 
sity profile returns to the same depth averaged value for 
the same surrounding conditions but with higher devia- 
tions at lower densities. The comparatively large errors 
relative to the number of the datapoints from the non- 
linear fit routine for the parameter Ue refiect a significant 
correlation between the 3 fit parameters (as could have 
been read off the covariance matrix) and hence judging 
the capableness of the applied model. 

From Fig. [iT] one can recognize once again the con- 
spicuous connection between the deuteron density and 
the screening energy like in the data of [^, [13] • High 
densities are Hnked to low screening energies because of 
moderately thick metal oxide layers as in the third case 
of [li. Fig. 13.(e,f)], Fig. [Blja). Examples are the ele- 
ments of the groups 3A (21 Sc, 39 Y and the lanthanoides 
Z = 57 - 71) and 4A (Ti, Zr and 72 Hf) emphasizing the 
chemical kinship with regard to the described surface re- 
actions in [l^ . Sec. 4.1], Sec. IIIII Low densities generate 
high screening energy findings due to shifts in the den- 
sity profile either in thick metal oxide layers or materials 
with low hydrogen binding ability as in [H, Fig. 13. (a- 
d)], Fig.[6ljb,c). Such can be recognized at the transition 
metals (groups 6A-8A: Z = 24 - 28, 42 - 46, 74 - 78) for 
example. It is argued in [l3] that the large enhancement 
findings are most likely due to electron screening because 
the data could be fitted well with the screening param- 
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eter f/fi. In view of the dispersion of the data points in 
[3, Fig. 1] their functional progression can also be 
described with the target model of [la, Sec. 4.3]. It imple- 
ments a simple static stepfunction for the density profile. 
The model can mimic a exponential like increase towards 
low energies quite passable by an inhomogeneous density 
profile with a super-deuterated surface layer alone with- 
out screening enhancement, i.e. Ue = 0, (igI . Fig. 10. c]. 
A existing screening increase can also be largely exag- 
gerated by the density profile of a deuterated zone in the 
metal with a limited thickness [B, Fig. lO.d]. Those were 
only static density profiles. A density profile dynamically 
changing with the energy as vindicated by Fig.[6l (igI . Fig. 
13] could perfectly imitate the exponential-like screening 
enhancement given the data distribution. In contradis- 
tinction thereto our data does not allow for such a de- 
scription as quantitatively demonstrated in (igI . Sec. 5, 
Fig. 12]. The Monte-Carlo code Srim for the simulation 
of ion stopping processes in matter was used in order to 
ratify the assumption of a homogeneous depth distribu- 
tion of the deuterons over the range of the ions [ll| . But 
Srim does not take into account the ability of hydrogen 
to diffuse. The homogeneity assumption was experimen- 
tally reconfirmed by a subsequent off-line ERDA on a 
4 MV tandem accelerator with the outcome that the dis- 
tribution is uniform within 10% for 'most' materials [l^ . 
Self-evidently a subsequent examination cannot detect 
dynamic changes but only the state of thermodynamic 
equilibration. Having been pointed to the problem of ox- 
idation [iHl RBS analysis was performed on the targets 
with the result that there were 'no detectable surface con- 
taminations' with the exception of Al where there was 
an AMDs layer with a thickness of about 150 monolayers 
0, Ell. Those findings prove that the resolution and 
sensitivity of the applied analysis techniques are too low; 
at least the passivation oxide layers from the unavoid- 
able exposition to air with the used equipment should 
have been visible. For both ERDA and RBS it is valid 
that light projectile ions with a kinetic energy of some 
McV cannot provide a wide energy spectrum of the ejec- 
tiles which would be necessary in order to resolve single 
atomic layers. Therefore a HIERDA with incident ener- 
gies of the heavy ions in the 0.1 GeV order of magnitude 
would be required with sophisticated magnetic analyz- 
ing systems (e.g. [Hi). This is additionally complicated 
by the circumstance that these methods deliver expres- 
sive results only if heterogeneous samples are made up of 
well defined layers. This is not fulfilled for the implan- 
tation targets with indistinct chemical composition and 
surfaces fractalized by embrittlement and beam deterio- 
ration (Fig. m [igI . Fig. 8]). So the applied methods are 
not able to detect metal oxides with a thickness of a few 
tens monolayers (some nanometers) , which is already suf- 
ficient to obliterate the screening enhancement (Sec. Illlj 
[l^ . Sec. 4.3]) while they are not thick enough to affect 
the applied density determination at 30 keV significantly. 

The thick oxide layer found on Al was defined to be of 
natural origin due to the said property of Al to readily 



oxidize on air. Hence a Kr ion sputtering treatment at 15 
or 35 keV was applied prior to the implantation measure- 
ments in order to remove those natural metal oxide layers 
which is the main difference from [H to [H and [li|. 
This procedure does not take into account that the ma- 
jor cause of the oxidation is contributed by the water in 
HV systems under deuteron irradiation which keeps going 
on nevertheless. While the high sputter yield of the Kr 
ions may allow for a surface cleaning the large Kr atoms 
thoroughly destroy the crystal structure of the target and 
get trapped in the material fractalizing the surface and 
thus possibly even promoting the oxidation process under 
subsequent deuteron irradiation since the necessary an- 
nealing step is omitted. The deviations in the screening 
energies between [l3|, [l3| and [l3| are in both direc- 
tions, anyhow giving an indication for the magnitude of 
the true error in the determination of the screening ener- 
gies in this way similar to our experiments on Ta (Fig.[5|). 
Whether the increase or decrease of the screening finding 
comes from an increase or decrease of the thickness of the 
oxide layer, or low hydrogen binding ability of the metal, 
or a too thin overheated target foil can scarcely be told 
afterwards on the basis of the available information. But 
there are peculiarities. It becomes not clear which beam 
currents were used, i.e. 54, 5 or 2.4 |^A, and how they in- 
fiuence the stability and the inferred screening values. In 
[l^ it was reported about instable yields dependent on 
the beam current for In and 'other elements with a low 
melting point' [8l|. The elements of the group IB (Cu, 
Ag, Au) had a small screening value in [ll.Tl3l| complying 
to the gas target value which became large in [l3|. That 
is in contradiction to the very low screening energy for 
Au of [§| and our finding of no screening. The extraordi- 
nary high screening value for Pd does not change with the 
Kr sputtering but matches best to the PdO value of [l3| 
which is another proof for the nevertheless continued ox- 
idation process. Due to the moderately thick metal oxide 
layer and the stable deuteron density close to the stoichio- 
metric ratio the metals of the groups 3 A, 4 A including 
the lanthanoides neither allow for a real screening obser- 
vation nor a simulated screening by density dynamics in 
the experiment of [l3|. In the coextensive publications 
[25l . [26| these metals where heated to 200° C thus over- 
coming the chemical bond between the metal atoms and 
deuterium conveying it into segregation and leading to 
a density drop of two orders of magnitude, i.e. the case 
of [H, Fig. 13.(e,f)], Fig.[6l[a) is transformed to [H, Fig. 
13.(a-d)], Fig. [6Ub,c). The than observed high screening 
energies again can be informally explained by the den- 
sity dynamics due to the high mobility of the deuterons 
induced by the high temperature and conjectural pro- 
moted metal oxide layer formation. This is made clear 
at the example of Ti where five datapoints taken at dif- 
ferent temperatures show the transition in [25l . [26l . Fig. 
3]. In [IJI it was reported on difficulties in attaining 
stable reaction yields for Ta at high temperatures which 
was subsequently not further elucidated. The stability 
test for the density in [2^1 is inapplicable since the used 
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analysis method cannot recognize the short time changes 
of the density. 

The intention of these comprising experiments was to 
find a connection between the observed screening energy 
and some electronic properties of the elements, some- 
thing that is to be underscored since it is an important 
step towards the understanding of this phenomenon. The 
authors propose the Hall coefficient to be this quantity 
stating that the effective density UeS of the free charge 
carriers, i.e. electrons and holes Hkewise, form a Debye 
sphere Rd around the deuterons and thus generate the 
screening potential (l3| : 



Table 11: Spearman rank correlation tests 
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with Pa the number density of the atoms T the temper- 
ature of the free electron gas and Zp the atomic number 
of the projectile. The classical Debye screening is, how- 
ever, not applicable for low temperatures (electron ener- 
gies below the Fermi energy) and dense plasmas (soHd 
states) where the quantum mechanical effects dominate 
and the screening effect depends only on the charged par- 
ticle density and not on the temperature [J, [H, [131 . Ad- 
ditionally, the motion of the bound electrons simulating 
the hole is not free but governed by quantum mechani- 
cal tunnehng between neighbor atoms. The fact that the 
screening energy is vanishing for high deuteron densities 
is explained by the assertion that these metal hydrides 
are insulators. This is not right for the majority of the 
metal hydrides which are metallicly or covalently bound 
and retain their metaUic properties. In fact, the elec- 
trons of the hydrogen are added to the conduction band 
of the metal. The Baranowsky-curve of the electric re- 
sistance of metal hydrides shows that the resistance at 
the chemical stoichiometric ratio is even lower than for 
somewhat lower densities and comparable to the metal 
[s^ l- Using a '^He beam on a deuterated rgPt target via 
the reaction d('^He,p)*He a screening energy was in- 
ferred about twice as high ((730 ± 60) eV at 1 — 3 |J.A) 
as for the d beam ((440 ± 50) eV) which was regarded as 
a confirmation of the Zp dependency l|19p of the Debye 
hypothesis In [l3| however the screening energies 

for ^He ((680 ± 60) eV) and d beams ((670 ± 50) eV) at 
Pt became equal without explanation. The inconsistency 
of the Pt-data also comprises the measurements [25;, |20| 
for the verification of the temperatur dependence (fTOj) of 
T~2 . With the exception of the room temperature data 
point the other four data points are equal within their 
error interval. So the temperature dependence is based 
on a single uncertain point. Furthermore, the findings for 
the metals of the groups 3A and 4A are in contradiction 
to it which cannot be resolved by the introduction of a 
highhanded function [H,!!!, Eq. (4)]. 
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, Table 1] and n^ff for the omitted elements. 
''Additional consideration of the sign of the Hall coefficient. 



In order to arrive at a more quantitative assessment of 
the hypotheses, methods of statistical data description 
and analysis can be applied (e.g. [i3|). Of the 58 exam- 
ined elements in [3, Table 1] the effective charge densi- 
ties calculated from the Hall coefficient Rh are selectively 
specified for the 25 elements with high screening values 
only, because the authors erroneously precondemned the 
others to be insulators with zero charge carriers. Ef- 
fective charge densities for the elements In, Sn, Sb, Pb, 
Bi which do not fit in the explanation scheme were also 
omitted fs^]. It needs to be particularly pointed out that 
the authors impute themselves a fiat error of 20% for 
the Hall coefficients but not the original experimental 
uncertainties. Both have severe impact on the interpre- 
tation. A visual survey based on Fig. [12] already shows 
that the distribution in the scatter plot for the case of 
Ue ^ n-cs is rather dispersed while the distribution for 
Ue ^ X roughly indicates a hyperboHc connection. The 
temperature dependent data points for Ti jH, [H, Table 
1] are additionally included in the bottom scatter plot 
(tagged with a star) demonstrating the transition from a 
high stable density to a low instable density allowing for 
the density dynamics which simulate the high screening 
findings. Three testing methods for continuous variables 
were used: Pearson's Hnear correlation r which assumes 
a Hnear association between the variables. The Spear- 
man rank correlation measures the monotone associa- 
tion between the variables and is therefore invariant un- 
der monotone transformations. Kendall's r is even more 
nonparametric since it uses only the relative ordering of 
the ranks by counting the inversions in the paired data 
points. It also enables the easy inclusion of errors by 
adaptive binning. The latter both are robust in opposi- 
tion to the linear correlation. All such tests attempt to 
falsify the null hypothesis of no correlation. Their corre- 
lation coefficient describes the strength of the correlation 
ranging in [—1, 1] where stands for no correlation and 
(-)l for total (anti)correlation. Complementary the P- 
value determines the significance of the obtained correla- 
tion, the lower P the higher is the significance. The re- 
sults are in tendency congruent wherefore representatives 
from the Spearman rank correlation are Hsted in Table [TTl 
The assumed functional dependency in l|19p is tested on 
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Table III: Regressions on Inrieff ^ In Ue 
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Figure 12: Scatter plots. Top: Ue ^ Jioff with UbH from the 
Hall coefficient Cnaii. Bottom: Ue ^ x with x being the 
deuterium contents in the compound MDa;. The star denotes 
temperature dependent points for Ti. 
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Sn, Sb, Pb, Bi is regarded (#4), with a concurrent de- 
crease in the significance. The consideration of the sign 
of the Hall coefficient lets the correlation approach zero 
(#5). On the other hand the enlarged data set has only 
sHght impact on the correlation to x (#6). Logarith- 
mizing (fT9| leads to a Hnear model with a slope of 6 = ^ 
and a positive intercept a containing the other quantities. 
Regression attempts based on it are listed in Table llllj 
(Ji is here the corresponding standard error of the two 
parameters. The infiuence of the error 5ncfi can be seen 
in comparison of #1 and #2, where #1 has the willfully 
by the authors imputed error of 20% while #2 adopts a 
more reahstic value of 10% in contrast. Anyway, the re- 
sulting slopes definitely stay behind the necessary value 
of i. Making worse the Q-value which is the goodness 
of the fit remains tiny. Conventionally if the goodness is 
smaller than 10"'^ the model is considered incorrect or 
the errors are still roughly underestimated. Here ( #3 ) 
the inclusion of the temperature dependent data [2 a . l26l . 
Table 1] and for the omitted elements In, Sn, Sb, 
Pb, Bi into the calculations leads to a slope close to zero 
and a goodness disqualifying the linear model. The value 
of the Hall coefficient of Pd was doubted in [l3| and re- 
placed by an own measurement with minor impact. So 
in both instances, correlation and regression, the expla- 
nation by the Debye hypothesis is ruled out. The dis- 
afSrmation of the Debye hypothesis on basis of the cor- 
relation tests alone might be disputable, together with 
the other aforementioned points it is decrepeted. The 
density hypothesis as an alternative clearly could not be 
falsified. Complemented by the preceding argumentation 
and the physico-chemical effects as described in [T^, Sec. 
4, 6]and Sec. [llllthe deuteron density dynamics provide 
the explanation for the alleged screening results. Thus 
the assumptions for reactions with heavier nuclei and ra- 
dioactive decay are refuted as well with the consequence 
that any experimental evidence offered for them need to 
come under scrutiny, Sec. IV BI IV CI 



the restricted data set (T3 . Table 1] (#1); the correla- 
tion coefficient is below 0.5 with a low significance. In 
contradistinction thereto the correlation to the density x 
has considerably higher values with utmost strong signif- 
icances (#2). Including the temperature dependent data 
[isl . HB, Table 1] into the calculations leads to a consid- 
erable decrease in the correlation coefficient (#3) which 
is further reduced when rics for the omitted elements In, 



B. Experiments with heavier nuclei 

As a consequence of Sec. IV Al the erratic high screen- 
ing findings of the other groups cannot really serve as a 
confirmation of our results. Therefore the experiments 
on Li target nuclei achieve special significance as a inde- 
pendent reassurance for the overall effect of the enhance- 
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ment of nuclear reactions in metallic environments. With 
regard to the reliability the most outstanding result dis- 
cussed in Sec. [V] is from the experiment d+Li in binary 
alloys with Pd and Au [l^. The a-yield from the re- 
actions at ^■'^Li (natural abundance: 7.42% ^Li, 92.58% 
''Li) in both alloys was observed normalized to the yield 
at 75keV in complete analogy to the procedure for the 
d+d reaction [l3| and set in relation to LiF targets. 
The inferred screening energies are (1500 ± 310) eV for 
PdLij: and (60 ± 150) eV for AuLi^:, where x was ini- 
tially at 5-10% for the alloys. The screening energy for 
PdLij: is therewith one order of magnitude higher than 
our value of (190 ± 50) eV for ^LiF which is in agreement 
with the simple theory The latter screening en ergy 
is substantially smaller than the (380 ±250) eV from [49l |. 
This is because a significant share of the increase of the 
observed S-factor towards low energies is caused by a 
2+ subthreshold resonance at E^, (*^Be) = 22.2 MeV. It 
needs to be included in the data analysis by sophisti- 
cated nuclear reaction theoretical calculations (isl . [50| . 
The advantage of this experiment on Li targets over deu- 
terium is that the small hydrogen atoms have a mobility 
in metals which is several orders of magnitude greater 
than other atom species. Thus, different to the auto im- 
planted deuterium targets there is no such fatal target 
atom density dynamics possible as in [l^, Fig. 13.(a-d)], 
Fig. [6l[b,c). The also here inevitable oxidation process 
will only decrease - but not completely supplant - the 
Li fraction in the surface layers. This explains together 
with the higher sputtering yield for light atoms the ob- 
served asymptotic bisection of the yield with the ion dose 
[igl . Fig. 2] at the monitor energy of 75 keV. Both leads 
to an inhomogeneous depth distribution of the Li target 
atoms with a lower Li fraction at the surface and due to 
the higher beam energies > 30 keV with lower impact. So 
this time the observed enhancement can be regarded as a 
lower hmit, too. Whereas the inferred value of the screen- 
ing energy needs to be corrected for the influence of the 
subthreshold resonance and the same questions regard- 
ing the screening energy calculation of the authors apply 
as in the case of deuterium. The low value for AuLi^; 
is so far in conformity with the negative findings for the 
d+d reaction in Au [13, [H, [l3|- A similar experi- 
ment [131 was later performed using the proton induced 
reactions on ^'^Li in an environment of Li2W04, Li-metal 
and PdLia; (x = 1%, 0.01%). The results for the screening 
energies of the reaction ''Li(p,a)Q; are (185 ± 150) eV for 
Li2W04, (1280±60) eV for the metal and (3790±330) cV 
for PdLii% which were obtained using standard proce- 
dures (dll), mi), [ii|,|4i|. The results for LiPdo.oi% and 
the reaction ^Li(p,Q!)'^He agree within Itr. In a micro- 
scopic view it is universally valid that the screening effect 
depends on the impact of the electronic configuration of 
the environment on the Coulomb barrier of the entrance 
channel only (e.g. (sil. Is^). i.e. the pure Coulomb energy 
is modified by a Yukawa factor for simplicity W{r) = 
_j_£p£te_^ with Aa being the screening length. As 
such the inferred screening energy is merely the second 



term in a Taylor-expansion of M^(r), i.e. Ue = 
and a coarse mathematical parametrization in the simple 
model |2|, [H, Eq. (20)]. The screening modifica- 
tion of the Coulomb potential only acts as if the projec- 
tile gained Ue- So there is no 'acceleration mechanism' 
in reality and one must neither decompose the screen- 
ing effect nor transfer the result of one environment to 
another as in [131 where the 'atomic' screening energy 
for ''Li— >H2 is used as a linear addend in the screening 
energy for p—*- (Li- metal or PdLia;). Consequently, the 
screening energy is independent of the isotopes in the 
reaction and should be equal for ^'^H+^'''Li in [l3| and 
[20| . Whereas there are two discrepancies between jlO] 
and [l3|: First, the more than twice as high screening 
energy for PdLi^; of Q relative to But 4 of the 

7 datapoints lay offside the fitted curve and only the fit 
error is given for Ue [20l . Fig. 1,2]. Second the assertion 
in [iOl that the yield remained stable better than 10% 
while [l3| observed a bisection of the yield at 6 C which 
is plausible due to irradiation effects. Both discrepan- 
cies can be explained by the different target fabrication 
techniques. In [l3| Pd and Li are made into an alloy by 
arc melting while in [lOl Li was inserted in a Pd-disk in 
a plasma discharge. The latter is prone to depth inho- 
mogeneities. This was verified by a NRA analysis of the 
target using the Ea = 958 keV resonance with a width 
of r = 4keV in the reaction ^Li(a,7)^^B yielding the 
ascertainment of a homogeneous depth distribution [2o| . 
However, the depth resolution of this method is limited 
by the energy uncertainty and spread of the beam and 
the width of the resonance. The most prominent exam- 
ple of the NRA is the i?N = 6.385 MeV resonance with 
a width of r = 1.8 keV in the iH(i5N,a7)i2C for the 
investigation of hydrogen distributions. It has a mini- 
mal resolution ranging 5 — 15 nm [H3|. So the resolution 
of the Li-NRA is worse given a 2.2 times higher width 
of the resonance. Since most of the yield is contributed 
by the topmost atomic layers here too ([l3, Sec. 4.3, Fig. 
10. d]), an enlarged Li contents below the NRA-resolution 
at the surface would explain the more than two times 
higher screening energy and the much lower decrease of 
the yield with the ion dose. The high screening value 
for PdLia: was regarded as a confirmation for the Debye 
model. If this was true the measurements for AuLi^; [19| 
should also have yielded a high value and not one close 
to zero, since the d+d screening energy for Au of 
was about 280 eV in difference to [3, 0, [H, and our 
observation. 

The theoretical model of the electron screening pre- 
sented in Sec. IIV Bl predicts different screening energies 
for different target material environments. In the case of 
an insulator the electron screening should reach the value 
of 190 cV, which results only from the gain of the elec- 
tron binding energies. For metallic environments the con- 
tribution coming from free electrons has to be included 
additionally. Due to different electron densities for Pd 
{rs = 1.4) and Li (rg = 3.4) the free electron contri- 
butions to the screening energy is equal to 660 cV and 
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420 eV, respectively. Thus, we finally expect total screen- 
ing energies of 190 eV for an insulating target material, 
610 eV for metallic Lithium targets and 850 eV for the 
PdLij, alloy. Experimental results, despite large uncer- 
tainties, confirm different electron screening energies for 
insulating and metallic materials with various electron 
densities. 

Extending this thread, a first effort was undertaken 
in [m to study the environmental influence for heavy 
nuclei using the (p,n) reaction on and ^''^Lu nuclei 
in an oxide, as pure metal, and as an alloy with Pd in 
the energy range 0.75 — 1.5 MeV. Because of insufficient 
cross section data the screening energies were obtained by 
comparison with the metal oxides VO2 and LU2O3. The 
inferred screening energies are (27 ± 9) keV and (33 ± 
11) keV for V and PdVio% and (32 ± 2) keV and (33 ± 
2)keV for Lu and PdLuio%. The comparison was done 
by taking the ratio of the yields between the metal and 
the oxide alias the insulator 

where is the efficiency of the neutron detector. It is 
now assumed that the ratio of the stopping cross sections 
between the two materials can be expressed by an energy 
independent constant a = ei{E)/eiri{E) which is mathe- 
matically doubtful considering Bragg's rule [HJl- So the 
following substitutions where done s„i{E) = a~^ei{E) 
and am{E) — f{E)ai{E) with the enhancement factor 
/ as in ([7]) with the presupposition of a constant S. 
The ratio of the integrals is further simplified by the en- 
ergy differentiation of the yields using the effective en- 
ergy as in dm), 111]) [nl, Eq. (5), (7), (8)] arriving at 
R{Ep) = af{Ecs). The screening energy resulted to- 
gether with a from a fit to the yield ratios. This proce- 
dure was, however, only applied to V. The screening ener- 
gies for Lu were gathered from the shift of the Lewis peak 
along the energy axis between the different targets orig- 
inating from a narrow resonance close to Ep = 0.8 MeV 
[21I . Fig. 3]. The Lewis effect comes from the discrete en- 
ergy loss of the projectiles in the target [5^. This energy 
shift was indeed erroneously interpreted as the screen- 
ing energy. As already pointed out the screening effect 
is merely a modification of the Coulomb barrier and no 
real energy shift. So this shift can not originate from 
the screening effect. Thus, it is probable that the energy 
shift is caused by target properties. Strikingly the oxide 
targets being the normalization standard are made by 
pressing a metal oxide powder into a cylindrical hole of 
a Cu disk. It is well known from powder metallurgy and 
silicate technology that pressing of a powder like in this 
case is insufficient in order to remove the hollow spaces 
between the powder particles unless a sintering step is 
performed. So the used metal oxide targets contain hol- 
low spaces with a size of the same order of magnitude as 
the powder particles. Consequently the stopping of MeV 
protons is heavily altered in comparison to a monolithic 



metal oxide ceramic and different in its mathematical de- 
scription to js^]. We observed effects of porous targets 
on the stopping [13, Fig. 2]. So the shift of the lewis 
peak can be explained by the differences of the stopping 
between the porous metal oxide target and the metal 
targets. Additionally, it is well known that the position 
and form of the Lewis peak depends very critically on 
the coraposition, homogeneity and contamination of the 
target [HsI . Is^ . IHtI also [i^. This in turn casts serious 
doubts on the results for V. A critical point of the data 
analysis ((20|) is the presupposition that the ratio of the 
stopping cross sections of the metal oxide and the metal 
is a constant over the energy. This is inappropriate for 
the porous target and can lead to a misinterpretation of 
the data. The conspicuously high errors of the screen- 
ing and a values from the fit — about 33% making the 
effect compatible with zero within 3cr — are a strong 
indication for a high correlation between the two fit pa- 
rameters showing the improperness of the fit model. The 
covariance matrix of the fit parameters could have given 
information about this. 

From the theoretical point of view the large screening 
energies obtained for the d+d reactions at energies below 
20 keV cannot be used for the estimation of the screening 
energies in the above case since the proton energy is much 
higher and does not fulfill the adiabatic approximation. 
Since the surrounding electrons are much slower than the 
protons, the resulting screening energy obeys rather as- 
sumptions of the sudden approximation and thus should 
be of order of a few keV in contradiction to statements 
involved in [23,[2i|. 



C. Radioactive decay of embedded nuclei 

As the electron screening enhances the cross section 
at low impact energies, a similar effect can be expected 
for the radioactive decay. However, since the energies of 
the decay products are fixed by the Q value, only a few 
nuclei with lowest-energy emitters are candidates for a 
measurable change in the lifetime. In general, for pos- 
itive charged ejectiles (a and decay), screening re- 
duces the Coulomb barrier and therefore enhances the 
decay rate while the opposite is true for /3~ decay. As re- 
cently pointed out by Zinner [HHl the effect of a changed 
Coulomb barrier is partially canceled by a modified Q 
value that stems from the extension of the screened po- 
tential into the inner part of the nucleus. For heavy nu- 
clei the effect can still be strong as the screening poten- 
tial scales approximately with the product of the charge 
number of the end nuclei. 

Recently, based on an extrapolation of the Debye- 
Hiickel electron screening model to low temperatures, it 
has been suggested that half-lives of radioactive isotopes 
may change by orders of magnitude if they are embedded 
in a metal lattice and cooled to cryogenic temperatures 
[2l|,[23,[2l,[23,[2i Hi. In support of these predictions, a 
series of measurements has been published the results of 
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Table IV: Decay of radionuclides embedded in host metals 



Ref. 


nuclide 


decay mode 


host 


prediction 


measurement 


[22] 


^^Na 


90% P+ 


Pd 


11% 


(1.2 ±0.2)% 


[59] 




100% (3- 


Au 


-34% 


(-4.0 ± 0.7)% 


[60] 


210p^ 


100% a 


Cu 


3300% 


(6.3 ± 1.4)% 



which are listed in the table IIV| together with the half- 
life changes predicted by the Debye-Hiickel model. The 
striking disagreement with the predictions have been at- 
tributed to an oxygen layer build-up on the metal surface 
leading to an insufBcient implantation of the radioisotope 

0- __ 

A recently published measurement |6ll.l62l| where ^^Na 
was activated in Al (and therefore deeply implanted) 
clearly shows a zero effect on a level of 0.04%, again 
in striking disagreement to the results by [l^l with a 
reported lifetime change of (1.2 ±0.2)% (see table lIVp . 
No description has been given in [131 how the data have 
been analyzed. If the 511-keV annihilation line has been 
included in the analysis, the results are certainly not cor- 
rect (see (s^l). For the a decay, even the observed 6% js^l 
change is surprising as embedding radioactive nuclei in 
metals and cooling the samples to cryogenic temperatures 
is a routine procedure in low temperature nuclear orien- 
tation (LTNQ) experiments since several decades [sH ]. 
Stone et al. [63| studied in detail the expected effect 
w.r.t. a decay on complete decay chains starting with 
^^''Rn, ^^^Ra, and ^^^Ac and compared it with available 
LTNO data. None of those data indicate any change of 
the lifetime of any of the nuclei involved when they are 
implanted into Iron, neither at room temperature nor 
when cooled to 20 millikelvin. The same applies for f3 
active nuclei in multiple host metals, see [63| and refer- 
ences therein. The precision of these measurements is 
typically 1% and less. Another follow-up measurement 
performed at ISOLDE/CERN focused on a possible 
change of the ^^^Fr {a decay) half-life when embedded in 
a metal and an insulator; there is also no clear effect (50% 
error) on a level of 0.3%. Severijns et al. [gBI investigated 
the a decay of ^^'^Es in Fe between 4 K and 50 mK and 
could not observe any effect on a level of 2%. Finally, also 
the P~ decay of ^^*Au embedded in Au and Al-Au has 
been measured independently by three different groups 
[HI, [13, [fill, and no lifetime change could be observed 
on a sub-percent level when the sample was cooled to 
«10 K. The latest result [g^I was measured with a 30 
times better error but the same conditions as in (59| . 

In conclusion, all the follow-up measurements are in 
agreement with the theoretical expectations presented 
already in Sec. IIV Bl The Debye-Hiickel screening can 
be applied only for temperatures higher than the Fermi 
temperature being typically 10^ K, far above the evap- 
oration temperature of metals. For lower temperatures 
one should not observe any temperature dependence of 
the screening energy. 

However, another effect can be expected a change 



of the lifetime by just embedding the unstable nuclei into 
a metal, but this requires an absolute measurement of 
the lifetime and therefore much more experimental ef- 
fort. A re-analysis of past lifetime measurements data 
with respect to the chemical composition could also re- 
veal such a dependence. For instance, the lifetime of ^^®U 
has been determined with electroplated samples (metaUic 
uranium), with UaOs, and other compounds, see (69l.[70|. 
Although the measurements scatter by 1-2%, no system- 
atic enhancement of the decay rate can be seen for the 
metallic uranium. As can be seen from the aforemen- 
tioned CERN measurement (6^ the effect will be small 
anyway even for high screening values but an evidence 
would be a great contribution to a better understand- 
ing of the screening mechanism from a very different ap- 
proach. 



VI. CONCLUSION 

We presented some new experimental electron screen- 
ing energies for d±d reactions taking place in different 
target material environments. We applied a differential 
data analysis method which gains the maximum infor- 
mation from the raw data. The method is indepen- 
dent of the unprecise stopping power coefficients and 
the actual absolute value of the deuteron number den- 
sity in the targets. It enables the on-line monitoring 
of the deuteron densities and the observation of short 
time deuteron density profile changes. Thus, it allows 
for the recognition and rejection of measurements with 
unwanted shifts in the density depth distribution pro- 
file. Therefore, it adequately considers the special situa- 
tion of potentially highly mobile hydrogen in soHd states 
where neither a homogeneous nor a stable density distri- 
bution can be presupposed any longer. The problem of 
the density dynamics is entangled with the effects from 
the actual target composition, i.e. the undesirable density 
profile changes occur in targets with low hydrogen bind- 
ing ability, like many of the transition metals, at elevated 
temperatures and heterogeneous targets with metal oxide 
or carbon layers or different (relatively) thin metal lay- 
ers. The formation of metal oxide layers is inevitable in 
common high vacuum systems used in experimental nu- 
clear physics while the other unpropitious environments 
were produced deliberately. Thorough investigation of 
the contamination layer formation showed their momen- 
tousness and assured together with the differential anal- 
ysis method that our screening energy values ranging be- 
tween 190 — 320 cV represent lower limits. In addition the 
alteration of the inferred screening energies due to layer 
formation under beam irradiation depends on many pa- 
rameters. Logically it makes no sense to measure larger 
portions of the periodic table since any observed mate- 
rial dependence results from differences in the chemical 
reactivity and related physico-chemical properties for the 
contamination layer formation, unless this problem is re- 
liably solved. Different to the other two groups our high 



20 



screening energy results were achieved at high densities in 
the proximity of the chemical stoichiometric ratio clearly 
without evidence for short time density profile shifts. 
Whereas the high screening results of the other groups 
were exclusively attained at low densities yielded from 
the customary analysis of the total yields of the mea- 
surements which is bhnd for the then happening density 
dynamics. The target diagnosis methods are unusable 
because of their too bad resolution and off-line applica- 
tion. So the inferred screening energies are conjecturally 
simulated by the density dynamics. Utilization of explo- 
rative statistics to the data sets including the temper- 
ature measurements sustains this explanation while on 
the other hand the Debye-Hiickel hypothesis is clearly 
falsified. It is Hkewise falsified from the theoretical side 
since calculations performed within an improved dielec- 
tric function theory predict only a weak material depen- 
dence of Ue on the valence electron density. The quan- 
titative scale of the phenomenon is not yet understood, 
since our analytical model still fails to describe the values 
by at least a factor of 2. So further unidentified effects 
play a role. Consequently any conclusion based on the al- 
leged material dependence of the inferred screening ener- 
gies is premature. For it the precise determination of the 
screening energies is demandable which is only feasible in 
an ultra high vacuum system with pressures well below 
lO^^'^ hPa, where only hydrogen and noble gases are in 
the residual gas, and equipped with in-situ target diag- 
nosis techniques. We performed the first measurements 
under UHV conditions, whose results confirm the previ- 
ous measurements and the framework of surface physics 



and chemical effects [7l] |. 

Nuclear reactions with heavier nuclei embedded in 
metallic environments gave evidence for an alike en- 
hanced screening effect. However, there are analogue 
problems. The results for deuterated metals with '^He 
projectiles are contradictory, most probably due to 
deuteron dynamics. The data for Li nuclei are partially 
conflicting between the Tohoku and the Bochum group, 
which used different target preparation techniques, and 
can be attributed to inhomogeneous densities and inad- 
equate diagnosis techniques as well. The results, how- 
ever, conflrm theoretical predictions based on the dielec- 
tric function theory concerning the free electron density 
of the target material. The screening energy data for the 
heavier nuclei V and Lu were obtained from a compar- 
ison between a metal and a metal oxide powder target 
ignoring the hollow spaces in the powder and its strong 
influences on beam stopping, thus disabling conclusions. 

As discussed, the predictions of the Debye-Hiickel hy- 
pothesis given by the Bochum group for the tempera- 
ture dependence of the radioactive decay of embedded 
nuclei could not be verifled by their own experiments; 
the measured values are orders of magnitude below their 
predictions. Moreover, their experimental results are in 
contradiction to all other experiments, in particular the 
LTNO measurements of the past 30 years. A material 
dependence is conceivable though a small effect. Oth- 
erwise it would have been already discovered given that 
nuclei of importance for nuclear technology have been 
investigated in multiple chemical compounds including 
pure metals for decades. 
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